A&A manuscript no. 

(will be inserted by hand later) 


ASTRONOMY 




Your thesaurus codes are: 

08( 08.14.1; 08.05.3; 02.04.1; 13.25.5) 


7Wi3 

ASTROPHYSICS 
1.2.2008 





Impact of medium effects on the cooling 
of non— superfluid and superfluid neutron 
stars 

-^^Ji, Schaab^, D. Voskresensky^'^, A. D. Sedrakian^, F. Weber^-^, and M. 
(K^ iWeigel^ 

On - 

'^^nstitute for Theoretical Physics, Ludwig-Maximilians Universitat, Theresienstrasse 
^III, D-80333 Munich, Germany 

Osesellschaft fiir Schwerionenforschung GSl, P. O. Box 110552, D-64220 Darmstadt, Ger- 
TTT3,ny 

Moscow Institute for Physics and Engineering, 115409 Moscow, Kashirskoe shosse 31, 



O^lax-Planck-Society, Research unit 'Theoretical Many-Particle Physics' at Rostock Uni- 
OO 

YuaJfsity, Universitatplatz 1, D-18051 Rostock, Germany 

^^vluclear Science Division, Lawrence Berkeley National Laboratory, MS: 70A-3307, Uni- 
^sity of California, Berkeley, CA 94720, USA 

5 ; 

T^ieived ; accepted 

I ' 

o ■ 

.^Qj^tract. Neutrino emission from the dense hadronic component in neutron stars 

il&>subject to strong modifications due to collective effects in the nuclear medium. 
• ^ 

implement new estimates of the neutrino emissivities of two processes operating 

.5—1 

iatfthe nuclear medium into numerical cooling simulations of neutron stars. The first 
process is the modified Urea process, for which the softening of the pion exchange 
mode and other polarization effects as well as the neutrino emission arising from the 
intermediate reaction states are taken into account. The second process concerns 
neutrino emission through superfluid pair breaking and formation processes. It is 
found that the medium effects on the emissivity of the modified Urea process result 
in a strong density dependence, which gives a smooth crossover from the standard 
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to the nonstandard cooling scenario for increasing star masses. For superfluid stars, 
the superfluid pair breaking and formation processes accelerate mildly both the 
standard and the nonstandard cooling scenario. This leads to a good agreement 
between the theoretical cooling tracks and the rather low temperatures observed 
for objects like PSRs 0833-45 (Vela), 0656+14, and 0630+18 (Geminga). The ro- 
bustness of our findings against variations in both the underlying equation of state 
of baryonic matter and the used fast cooling processes is demonstrated. Hence we 
conclude that the two recalculated neutrino emissivities studied here enable one to 
reproduce theoretically most of the observed pulsar temperatures by varying the 
masses of neutron star models. 
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1. Introduction 

The physics of neutron star cooling is based on a number of ingredients, among 
which the neutrino emissivity of the high-density hadronic matter in the star's core 
plays an important role. Depending on the dominant neutrino-emission process 
at the early stages of the thermal evolution, the cooling simulations follow either 
the slow (standard) or the fast (nonstandard) scenario of thermal evolution. In 
the first case the dominant neutrino-radiation reactions are the modified Urea and 
the bremsstrahlung processes. In the second case, these are the pion (kaon) /?- 
decay processes, direct Urea on nucleons and hyperons, as well as their analogous 
reactions taking place in the deconfined quark phase. The main difference in the 
cooling efficiency driven by these processes lies in the rather different phase spaces 
associated with these reactions. In the first case the available phase space is that of 
a two-baryon scattering process, while in the second case it is that of a one-baryon 
decay process. 

Numerical simulations of neutron star cooling, incorporating these types of neu- 
trino emission, have been extensively performed in the past (see, for instance, Tsu- 
ruta 1966, Richardson et al. 1982, Van Riper 1991, and Schaab et al. 1996a). These 
calculations suggest that the combined soft X-ray data from Einstein, EXOSAT, 
and ROSAT are roughly consistent with the slow cooling scenario, depending on 
the equation of state and the uncertainties associated with the behavior of super- 
dense matter. Some pulsars, however, as for example Vela and Geminga, possess 
rather low temperatures and thus seem to call for a more rapid cooling than is 



obtained for the modified Urea process. The fast coohng scenarios, on the other 
hand, generally have the tendency to underestimate the surface temperatures of 
these pulsars. So a natural question, motivated by these observations, is whether 
the physics of neutrino emission from the star's dense core will be able to resolve 
this problem. This will be discussed in this paper. Obviously, the corresponding 
modifications must be such that the theoretical cooling curves describe both the 
hotter as well as cooler classes of pulsars, which could be linked to variations in 
gross-structure parameters, like the star's mass (Voskresensky & Senatorov 1984 
(VS84), 1986 (VS86)). 

The aim of the present work is to implement a number of medium-modified neu- 
trino emissivities, which dominate a neutron star's thermal evolution, in numerical 
cooling simulations of such objects. We shall demonstrate that by means of this 
one is indeed able to achieve agreement with both the high as well as low observed 
pulsar temperatures. 

It is a well established fact that neutrino emission at the early stages of the stan- 
dard evolution is dominated by the modified Urea and bremsstrahlung processes. 
The neutrino radiation in the modified Urea process was first estimated by Bah- 
call & Wolf (1965) and by Tsuruta & Cameron (1965). Bremsstrahlung processes 
dealing with the neutral currents were first evaluated by Flowers et al. (1975). 
Then the corresponding rates have been more closely studied by Friman & Maxwell 
(1979, FM79) resulting in sufficiently higher rates. The calculation were performed in 
the free one-pion-exchange approximation to the long-range part of the nucleon- 
nucleon (NN) interaction, supplemented by a parametrization of the short-range 
part of the NN interaction by the Landau Fermi-liquid parameters. Medium effects 
enter these rates mainly through the effective mass of the nucleons. Therefore the 
density dependence of the rates is rather weak and the neutrino radiation from a 
neutron star depends only very weakly on its mass. This is the reason why the stan- 
dard scenario based on the FM79 result, though complying well with a few slowly 
cooling pulsars, fails to explain the data of the more rapidly cooling ones. 

Here, we shall carry out detailed simulations of the cooling of neutron star mod- 
els, which incorporate the softening of the one pion exchange mode, other medium 
polarization effects, like the inclusion of the nucleon-nuclcon correlations in the 
vertices, as well as the possibility of neutrino emission from intermediate reac- 
tion states^ (VS84, VS86). We also include the processes of neutrino pair radiation 

^ E.g., the pion participating in the exchange between the nucleons may also radiate 
the neutrinos or decay in intermediate reaction states to the nucleon-nucleon hole with 



from superfluid nuclcon pair breaking and formation mechanism first estimated by 
Flowers et al. (1976), and then more closely studied by Voskresensky & Senatorov 
(1987) and Senatorov & Voskresensky (1987) using the closed diagram technique. 
The medium modifications of these rates result in a pronounced density dependence, 
which, in turn, links the cooling behavior of a neutron star decisively to its mass. 

Using a collection of modern equations of state for nuclear matter in its ground 
state, which covers both relativistic as well as non-relativistic models, we shall 
demonstrate the robustness of the new cooling mechanisms against variations in 
the equation of state of super-dense neutron star matter. Part of these variations 
are caused by the possible superfluid behavior of neutron star matter. To demon- 
strate this effect on the new cooling mechanisms, we proceed in three successive 
steps, starting from stars made up of non-superfluid matter. These models are then 
supplemented with superfluidity, with superfluid gap parameters taken from the 
literature. Unfortunately, the gap energies are not very well known because of un- 
certainties in the effective baryon masses, the nucleon-nucleon interaction, and po- 
larization effects (cf. Chen et al. 1993). To account for these uncertainties, their 
values are finally varied about the theoretically determined ones. Since there are 
also uncertainties concerning the possibility of fast cooling processes the calcula- 
tions are done with both the nucleon direct Urea and the /3-decay processes in pion 
condensates (pion Urea) and additionally without any fast cooling process. 

We shall show that the above mentioned medium-modifications lead to a more 
rapid cooling than obtained for the standard scenario. Hence they provide a possible 
explanation for the observed deviations of some of the pulsar temperatures from 
standard cooling. Particularly, it is shown that they provide a smooth transition 
from standard to nonstandard cooling for increasing central star densities, i.e., star 
masses. 

The paper is organized as follows. In Sect. 2 we briefly discuss the incorporated 
neutrino emission processes which comprise the modified Urea, superfiuid nucleon- 
pair breaking and formation, direct Urea and pion Urea processes. In Sect. 3 the 
neutron star cooling simulations and the physical input quantities used in the calcu- 
lations are described. In Sect. 4 we present the results of the numerical calculations 
and compare them with the observed data. The conclusions are summarized in the 
last section. 

subsequent neutrino radiation. This essentially modifies the absolute value as well as the 
density dependence of the modified Urea process rate. 



2. Neutrino Emissivities 

2.1. Medium Effects on the Modified Urea Process 

The modified Urea (MU^) process is linked to neutrino emission via a weak de- 
cay of a nucleon in the presence of a bystander nuclcon. The latter ensures that 
both momentum and energy are conserved in the reaction. The two topologically 
different diagrams contributing to these processes are schematically shown in Fig. 
1. The first diagram shows the neutrino radiation from a nuclcon "leg", i.e., from 
an initial or final -state nucleon. The shaded block represents the total nucleon - 
nucleon interaction amplitude in the medium (including the softening of the pion 
mode, corrections of the vertices by the nuclcon-nucleon-hole and delta-nucleon- 
hole loops, and modifications of the short-range interaction by nucleon- nucleon-hole 
and delta-nucleon-hole loops), while the filled circle stands for the nucleon-nucleon 
correlations. The second diagram illustrates the neutrino radiation from the interme- 
diate scattering states (like neutrino radiation by the pion and the nucleon-nucleon 
hole loop being present in intermediate reaction states). Such medium effects as well 
as the second reaction channel were not incorporated in the calculations by FM79 
and shall lead to a rather new cooling behavior, as we will see below. 



Table 1. Key to abbreviations 



abbreviation 


meaning 


MU-FM79 


modified Urea process calculated by FM79 


MU-VS86 


modified Urea process calculated by VS86 


NPBF 


neutron pair-breaking and formation 


PPBF 


proton pair-breaking and formation 


DU 


direct Urea process 


PU 


pion Urea process 



The main contribution to the nucleon-nucleon interaction amplitude at densities 
n ^ 0.5— 0.7 no, where no — 0.16 fm~^ is the nuclear matter saturation density, is de- 
termined by the in-medium one-pion-exchange supplemented by medium-modified 
^ We shall abbreviate the modified Urea process calculated via the free one-pion exchange 
model by Priman and Maxwell with MU-FM79. The abbreviation MU-VS86 will refer 
to the process which accounts for the additional included medium modification effects. 
The abbreviations of the various neutrino-emitting processes considered in this paper are 
summarized in Table 1. 





Fig. 1. Neutrino emission from a nucleon leg (left graph) and from intermediate scattering 
states (right graph). Details are given in the text. 



vertices (VS86), whereas the peculiarities of a more local part of the interaction play 
an important part only at rather small densities, n ^ 0.5no (Blaschke et al. 1995). 
At n ^ no the emissivity arising from the modified Urea is dominated by the sec- 
ond diagram of Fig. 1, whereas the first diagram gives only small corrections (VS86, 
Voskresensky et al. 1987, Haubold et al. 1988). The contributions of the first di- 
agram, which gives the only contribution to emissivity calculated by FM79, can 
therefore be neglected. The resulting emissivity reads in a simplified notation as 
follows (VS86, Migdal et al. 1990), 



^(MU-VS86) 



2.4 X 10^"* Fi I — 
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where Tg = T/IO^K is the temperature, and mjv are the effective and 
bare nucleon mass, respectively, PFn{n) is the density dependent neutron Fermi- 
momentum, and the factor 



4r2 Vno 



2/3 



is the correction due to the pion decay from intermediate states. It reduces to unity 
when the intermediate state pion decay processes are ignored and only contributions 
from the weak current decay in the nucleon nucleon-hole loop being present in 
intermediate reaction states are taken into account. 

The quantity F accounts for the nucleon-nucleon correlations in the ttNN ver- 
tices. The effective pion gap, d)^, expressed in units of m^r = 140 MeV, is given by 
^^Fn{n)\ ~ —D~^ (w ~ /i^, k = PFnin)), where is the in-medium pion Green 



function, /i,r is the pion chemical potential. (For simplicity, we set the chemical 
potentials of 7r+,7r~, and 7r° mesons equal to /i,r — 0.) 

The values of HP' and F for isospin symmetric nuclear matter at saturation den- 
sity are extracted from the atomic nuclei data. For the case of interest, /i^ — 0, they 
read F(no) 0.4 and ixl^ino) ~ 0.8-0.95 (Migdal et al. 1990). For highly asymmet- 
ric nuclear matter, especially for n > no, the magnitude and density dependence of 
uP'in) and r(n) are not well known. In the present calculations we approximate the 
form- factor in the effective ttA^'A'' vertex, r(n), by the formula 

« ] , u/3 ' (2) 

1 + 1.4 (n/no)'/^ 

which shows the explicit functional dependence on the density.^ The effective 
pion gap Co^{n) has been estimated in Migdal (1978), Voskresensky & Mishustin 
(1978, 1982), and Migdal et al. (1990) for different parameter choices. We apply 
two different parametrizations in the cooling calculations (see Fig. 2). The first one 
assumes a phase transition into the pion condensate phase at = 3no (Brown & 
Weise 1976, Migdal 1978, Migdal et al. 1990), while the second one assumes that 
there is no phase transition (for a more detailed discussion see Sect. 2.3). Due to the 
pion-fluctuation effect the pion condensation sets in via a first order phase tran- 
sition (see Dyugaev 1975, Voskresensky & Mishustin 1978, 1982). This manifests 
itself in a jump of the pion gap at n = from a positive value to a negative one, see 
solid line in Fig 2.^ It has to be emphasized that the density dependence of a)^(n) is 
rather unknown. Our choices for the parametrization are motivated by microscopic 
many-body calculations. Somewhat different choices for the values of a)^(n), r(n) 
and the critical density for the onset of pion condensation affect of course the result- 
ing cooling rates. The qualitative conclusions, however, remain unchanged. Indeed, 
as we will see below, the cooling curves for the two different parametrizations of 
(D^(n) do not differ very much (see e.g. Fig. 8). 

The value of the parameter a is given by a = 1 for n < n^, and a = \f2 for n > 
nj. This parameter accounts for the changes caused by the neutrino emissivity of the 
modified Urea process below and above the critical point where pion condensation 
occurs (cf. VS86). 



^ This implies in our case that the Landau-Migdal parameter g' scales as (n/no)^^^ 
leading to an increasing repulsion for higher densities. 

^ Without this behavior one would get an anomalous increase of the neutrino emissivity 
via bremsstrahlung in the vicinity of the critical point. This effect is analogous to the 
critical opalescence phenomenon (VS86). 
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Fig. 2. Effective pion gap Cj^ /m^ versus density with (solid line) and without pion con- 
densation (dashed line), no = 0.16 fm~^ denotes the nuclear saturation density. 

The factor 



takes roughly into account the suppression caused by the nucleon pairing, where 
A„ and Ten are the superfluid gap and the critical temperature for superfluid phase 
transition of neutrons. The factor C('^p) is defined in the analogous way and takes 
into account the proton-proton pairing. 

2.2. Cooper Pair-Breaking and Pair-Formation Processes 

When a neutron star cools down to temperatures in the vicinity of the critical 
temperature for pairing of neutrons and protons, pairing correlations play an in- 
creasingly important role in the dynamics of the star's thermal evolution. Generally, 
the onset of superfluidity tends to slow down the cooling rate of a neutron star, since 
the neutrino emission processes are drastically suppressed. It enables however two 
additional processes as noticed in Flowers et al. (1976) and Voskresensky & Sena- 
torov (1987): the superfluid neutron-pair breaking and formation process (NPBF) 
and the superfluid proton-pair breaking and formation process (PPBF). 

The superfluid in a neutron star can be considered as a two-component sys- 
tem, which, for a fixed density and temperature, consists of paired quasiparticles 
in the condensate and elementary excitations above the condensate. Their associ- 
ated quasi-equilibrium densities are controlled by Cooper-pair-formation and pair- 
breaking processes. These processes, which are rather frequent at temperatures in 




(3) 



the vicinity of Tc, become successively suppressed at lower temperatures, because 
of an exponential increase in the number of paired particles. As shown in Fig. 
3, these processes proceed with the emission of neutrino pairs via the reactions 
{NN} ^N + N + u + P and N + N ^ {NN} + u + P, where {A^A^} denotes the 
Cooper pair, N an excitation. 



Fig. 3. Neutrino emission from Cooper pair-breaking (left graph) and pair-formation 
processes (right graph). Details are given in the text. 

The respective emission rates can be calculated either in terms of the closed di- 
agram formalism, using the normal and anomalous Green's functions for the paired 
nucleons, or in terms of the method of Bogoliubov's transformations applied to re- 
action rates in the normal ground state. As far as one works in the quasiparticle 
limit, which is an excellent approximation for neutron star matter, the preference 
for one or the other method is primarily a matter of taste. However, in a broader 
context, the Green's functions method provides a more general treatment by allow- 
ing for off-mass-shell effects and a systematic diagrammatical representation of the 
approximations to the self-energy functions (scattering rates). A first estimate of 
the reaction rates was given by Flowers et al. (1976) using the Bogoliubov's trans- 
formations. It was recalculated by Voskresensky & Senatorov (1987) using closed 
diagram technique. Voskresensky & Senatorov have also included the contribution 
of the axial vector coupling and used different parameterizations of the weak cou- 
pling vertices. This results in a rate being higher by an order of magnitude. Besides 
the NPBF process, the PPBF process was considered, too. One finds in lowest order 
of the expansion parameters T/cpn and T/epp {cpn and epn are the Fermi-energies 
of neutrons and protons, respectively) for the neutron and proton components (see 
Voskresensky & Senatorov 1987) 
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It can be seen that in the hmit A/T ^ 1, the rate of these processes is ex- 
ponentially suppressed, as it is the case for the two-nucleon process considered 
above. However, because of mild phase space restrictions (we deal here effectively 
with a one-nucleon phase-space volume) these processes may considerably con- 
tribute to the neutrino emissivity of the neutron star. The maximum value of 
(A(T)/A(T = 0))^/(A(r)/r) depends only on the ratio A(T = Q)/Tc. For ^So- 
pairing the maximum value is equal to « 10~^, for "^P2-pairing it is w 10^^. This 
causes the emissivity of the NPBF process in the core (■^P2-paring) to be smaller 
by three orders of magnitude (for equal gap energies) compared to the emissivity 
in the crust (^So-pairing). 



2.3. Direct Urea Process and Pion Urea Process 

Along with the two processes described above we include in the numerical calcula- 
tions other relevant channels of neutrino emission, too. Though we shall see that the 
processes described above lead to "intermediate" cooling rates, it still appears use- 
ful to distinguish between "standard cooling" and "nonstandard' or "fast cooling" , 
depending on whether the predominant process is the modified Urea process or a 
faster one, as nucleon or hyperon direct Urea, processes in pion- or kaon-condensed 
matter, or those caused by the presence of quarks in dense matter. 

As examples of a fast cooling processes, we shall use the nucleon direct Urea 
(DU) and the pion Urea (PU) processes. These processes are rather representative 
for other fast processes. Indeed, all the fast cooling processes lead qualitatively to 
the same temperature dependence of the emissivity, T^. Moreover the values of 
the critical densities beyond which these processes are possible are close to each 
other. 

The DU processes 

n — > p + e~ -|- i^e and p + e~ — > n + i/g (7) 



can only occur if the proton fraction exceeds some critical value (about 11-13 % 
depending on the composition) in order to fulfill energy and momentum conservation 
(see Boguta 1981 and Lattimer et al. 1991). In contrast to nonrelativistic EOS's 
(e.g. UV14+UVII) which predict proton fractions that he below these critical values, 
relativistic EOS's contain sufficiently large proton fractions such that the direct Urea 
process becomes possible. However, the proton fraction depends crucially on the 
symmetry energy which is unfortunately not well known at higher densities. We shall 
use the emissivity calculated by Lattimer et al. (1991), however corrected by the 
suppression factor (see eq. (2)) accounting for the nucleon-nucleon correlations 
in the weak interaction vortex (for simplicity, wo assume the same suppression factor 
in all vertices, cf. discussion in Voskresensky & Senatorov 1987): 

e(DU) ^ 4 ^ 1027 M (!!}k\ " TiCiAn) C(A,) . (8) 

\no J \mN J cm'^ sec 

The impact of the direct Urea process on the cooling will be shown in Sect. 4.4, 
where the final results are discussed, whereas it is neglected in the other cooling 
curves in order to demonstrate the efficiency of MU-VS86 and NPBF and PPBF- 
proccsscs. 

Even more uncertain is the possibility of the hyperon direct Urea (see Prakash et 
al. 1992). Whereas the contribution of the hyperon direct Urea is small compared to 
the contribution of the nucleon direct Urea in the non superfluid case, the hyperon 
direct Urea might become important if the nucleons form superfluid pairs. If however 
the hyperons become also superfluid, then their contribution is again neglible. The 
superfluid phase transition of hyperons has unfortunately not been studied so far 
since it is rather difficult to implement superfluidity in relativistic treatments with 
hyperons. We therefore neglect this contribution to the neutrino emissivity. 

The pion condensation was suggested by Migdal (1971), Scalapino (1972), and 
Sawyer (1972) and then considered by many authors (e.g. Brown & Weise 1976, 
Migdal 1978, and Migdal et al. 1990). In the last years, some arguments against 
the occurrence of pion condensation in neutron stars were given in the literature 
(see, for instance. Brown et al. 1995). These arguments are generally based on 
the rather strong increase of the Landau-Migdal parameter g' with the nucleon 
density. However other effects which might be important for the pion condensation 
problem (cf. Migdal et al. 1990) were not incorporated. Thus one may conclude 
that the question whether a pion condensate occurs in neutron stars is still not 
settled (see Ericson & Weise 1988 and Kunihiro et al. 1993 for extensive reviews). 



We shall therefore use models both with and without pion condensates applying 
two parametrizations ofuj^{n) (see Sect. 2.1). 

Additionally the neutrino emissivity in kaon condensation, which is now favored 
by various authors (e.g. Brown et al. 1988 and Thorsson et al. 1995) as well as 
the critical density for the onset of condensation are quite similar to the case of 
pion condensation. The additional tan(8c)-factor in the kaon condensation case, 
where &c — 0.223 is the Cabibbo-angle, is partly compensated by a probably larger 
value of the suppression-factor due to a possible suppression of the correlations 
between strange and non-strange baryons. The cooling curves for the models with 
kaon condensation look quite similar (see, for example. Umeda et al. 1994 and 
Schaab et al. 1996a). We shall therefore study only pion condensation. 

Above threshold density for pion condensation n^, the contribution of eq. (1) to 
the neutrino emissivity of the MU process is to be supplemented by the correspond- 
ing pion Urea processes (cf. Maxwell et al. 1977, VS84). For the latter processes we 
use a simplified expression including the nucleon-nucleon correlation effect in the 
ttNN vertices (VS84, Migdal et al. 1990), 



for densities n > nj. Here the neutron Fermi momentum is expressed in units of 



developed condensate implies 6 — > 7r/2. We shall use an ansatz that interpolates 
between these two extremes. Since 7r+'~ condensation probably reduces the energy 
gaps of the superfluid states by an order of magnitude (see Takatsuka & Tamagaki 
1980), we assume, for the sake of simplicity, that superfluidity vanishes above ric- 
Finally we note that though the PU processes have genuinely one-nucleon phase- 
space volumes, their contribution to the resulting emissivity is suppressed relative to 
the direct Urea by an additional vertex factor due to existence of the additional 
ttNN vertex in the former case. 

The modifications of energy density and pressure due to the pion condensation 
are taken into account as 




X 



erg 



(9) 



m^, and sin^ ~ ^/2\^P\/rn^ for ^ <C 1, whereas at large densities (n ) a well- 




(10) 



where, again, lj'^ is given in the units m^. 



3. Stellar Composition and Structure 



The general relativistic equations of stellar structure and thermal evolution (cf. 
Thorne 1977) were solved via a numerical code based on an implicit finite difference 
scheme handled by a Newton- Raphson algorithm. Details can be found in Schaab 
et al. (1996a). The physical input quantities, are summarized in Table 2. Besides 
the new neutrino emission processes discussed in this paper, we include in the 
simulations the traditional neutrino processes too, which are discussed in greater 
detail by Schaab et al. (1996a). ^ 



Table 2. Input quantities used for the cooling simulations 



Parameter 


References 


Equations of state: 




crust 


Baym ct al. (1971), Negele & Vautherin (1973) 


core 


sec Tabic 3 


Superfluidity 


see Table 4 


Heat capacity 


Van Riper (1991), Shapiro & Teukolsky (1983) 


Thermal conductivities: 




crust 


Itoh et al. (1984), Itoh (1983), Mitake et al. (1984) 


core 


Gnedin & Yakovlev (1995) 


Neutrino emissivities: 




pair-, photon-, plasma-processes 


Itoh et al. (1989) 


bremsstrahlung in the crust 


Itoh & Kohyama (1983), Pethick & Thorsson (1994) 


bremsstrahlung in the core 


FM79, Maxwell (1979) 


MU-FM79 


FM79 


or instead: MU-VS86 


eq. (1), VS86 


NPBF and PPBF processes 


oq. (4,5), Voskrcscnsky & Senatorov (1987) 


DU-process 


eq. (8), Lattimer et al. 1991 


PU-process 


eq. (9), VS84 


photosphere 


Van Riper (1988) 



® The two-nucleon neutral current processes (bremsstrahlung-processes) are also modi- 
fled by in- medium effects (VS86). Their contribution is however much smaller than the 
corresponding contribution of the MU-VS86 process. These modifications are therefore not 
implemented here. 



3.1. Equation of State 

For the outer and inner crust we adopt the equations of state of Baym et al. (1971) 
and Negele & Vautherin (1973). The transition density between the ionic crust 
and the core of a neutron star is assumed to be ptr = 1-7 x 10^^ g cm~^ (Pethick 
et al. 1995). The high density matter in the cores of neutron stars is described 
by a collection of modern equations of state, which consists of a non-relativistic 
Schrodinger-based model as well as three relativistic field theoretical ones. The 
details of these models are summarized in Table 3. One of the significant differences 
between these equations of state is that the non-relativistic model treats neutron star 
matter as being composed of neutrons and protons only (which are in /3-equilibrium 
with leptons). The relativistic models account for all hyperon states that become 
populated in the cores of neutron stars as predicted by theory (Glendenning 1985). 
Dynamical two-nucleon correlations calculated from the relativistic scattering T- 
matrix in matter are contained in the relativistic Brueckner-Hartree-Fock (RBHF) 
equation of state. The underlying one-boson-exchange interaction is Brockmann's 
potential "B". At densities larger than about three times normal nuclear matter 
density it has been joined with a relativistic Hartree-Fock (RHF) equation of state 
(Huber et al. 1994, 1996). The non-relativistic treatment of the EOS leads generally 
to a softer EOS than the relativistic treatment. However, since in the non-relativistic 
case one takes only nucleons and leptons into account the relativistic EOS's with 
hyperons become almost as soft as the non-relativistic EOS's. Both treatments 
therefore mainly differ in the resulting composition. This can be important for the 
nucleon and hyperon direct Urea (see Sect. 2.3). The equations of state are shown 
in Fig. 4 and the corresponding neutron star sequences in Fig. 5. 

Table 

3. Collection of nuclear equations of state (EOSs). Abbreviations: nvar=nonrelativistic 
variational, RH=relativistic Hartree, RHF=relativistic Hartree-Fock, RBHF=relativistic 
Brueckner-Hartree-Fock. 

EOS Composition Interaction Method Ref. 

UVi4+UVn P,n,e,^ Urbana Vm+VII nvar Wiringa ot al. (1988) 

HV p,n,A, E, H,e,/i a,uj,p RH Weber & Weigel (1989) 

RBHF(B)-|-HFV p,n,A, S, H, A,e,/i a, ui, n, p, r), 5 RBHF, RHF Huber et al. (1994) 

G300 P,ii,A, E, H,e,/U <y^ijj,p RH Glendenning (1989) 
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Fig. 4. Equations of state (pressure versus energy density) used in this paper. The dots 
refer to the maximal central densities reached in the most massive neutron stars con- 
structed for these EOS's. The modifications caused by a pion condensate alter the EOS's 
only insignificantly on this scale and are thus not shown. 
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Fig. 5. Star masses in units of solar mass versus central star density for the EOS's used 
in this paper. 



3.2. Superfluidity 

The following superfluid regions inside neutron stars are taken into account: 

• 7 X 10^^ g cm~^ < p < 2 X 10-^^ g cm~^: neutrons in ^So-pair state, 

• 2 X 10^^ g cm~^ < p < 4 X 10^^ g cm~^: protons in ^So-pair state, and 

• 2 X 10^^ g cm~^ < p < 5 X 10^^ g cm~^: neutrons in ^P2-pair state. 



The proton gap is taken from Wambach ct al. (1991), those of ^Sq- and ^P2-paired 
superfluid neutrons are from Ainsworth et al. (1989) and Amundsen & 0stgaard 
(1985), respectively. Details are given in table 4. It should be emphasized that the 
extension of the superfluid density regime depends on the equation of state. The 
above values are computed for the HV equation of state, and Fig. 6 displays the 
profiles of the superfluid gaps for this particular model. 
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Fig. 6. Superfluid energy gaps computed for the HV EOS. We neglect superfluidity in the 
models where pion condensation is allowed for densities n > = 3no (dotted line). 



Table 4. Maximal gap energies, A^*™, critical temperatures, T™*", and density ranges, 
p, in neutron star matter computed for the HV equation of state 



Neutron ^So 



Neutron 



Proton ^So 



Reference 



0.62 
0.8 

" 5 X 10^5 



0.25 
1.6 



A,"^"" [MeV] 1.13 

p[gcm"^] 7x10" -2x10" 2x10 _ _ _ 

Ainsworth et al. (1989) Amundsen & 0stgaard (1985) Wambach et al. (1991) 



2 X 10" - 4 X 10 



,14 



Because the pairing gaps are not too accurately known, we shall also study 
cooling scenarios of superfluid neutron stars where the magnitudes of the gaps 
are varied in an ad hoc fashion about their theoretically determined values. This 
procedure may serve to reveal the dependence of our theoretical cooling scenarios 
on the uncertainties associated with the superfluid gaps. 



As already mentioned in Sect. 2.1, the superfluid gaps are reduced by an order 
of magnitude in tt"*"'" condensed matter (see Takatsuka & Tamagaki 1980). We 
therefore neglect superfluidity in pion condensed matter completely. In passing we 
emphasize that in a more elaborate treatment, one needs to pay attention to the 
effects of interactions between tt~^'~ and tt" condensates, simultaneous appearance 
of K~ and K° condensates (Brown et al. 1994, Kolomeitsev et al. 1995) which 
affect the neutron and proton gaps differently. Furthermore, neutron-proton pairing 
in the ^£'2-triplet state for almost symmetric nuclear matter^ may dominate the 
traditionally considered pairing (Aim et al. 1996). 

Finally when performing numerical cooling calculations for superfluid neu- 
tron stars, we follow the traditional procedure by suppressing the two-nucleon 
reaction neutrino emissivities, the thermal conductivity, and the heat capacity 
of the nucleonic constituents of the stellar matter by factors which behave like 
exp(— An (pj/A-bT) for T <C Tc, where A„(p) denotes the gap, T the temperature 
(see MaxweU 1979, Gnedin & Yakovlev 1995)^. 

4. Results and Discussion 
4-1- Observational Data 

Compact X-ray sources have been detected by the X-ray observatories Einstein, 
EXOSAT and ROSAT during the last two decades. Among them, 14 X-ray sources 
- observed at least by one of these satellites - were identified as radio-pulsars (with 
the exception of Geminga, which is known to be radio- quiet). The information 
obtained from these detections is not always sufficient to extract the effective surface 
temperature of the corresponding pulsar. Therefore this sample of pulsars is divided 
into three different categories (see Ogelman 1995): 

1. The detection of three pulsars (PSR's 1706-44, 1823-13, and 2334+61) contain 

too few photons for spectral fits. The specified luminosities are calculated by 

using the totally detected photon flux. These pulsars are marked with triangles 

in the figures that will be discussed below. 

^ In neutron star matter with hyperons based on relativistic EOS's the neutron and 

proton fractions tend to each over with increasing density and are almost equal in the high 

density limit (cf. Weber & Weigel 1989 and Glendenning 1989). One must distinguish this 

case from the nonrelativistic treatment without hyperons, where in the presence of kaon 

condensation isospin symmetric matter is preferred (cf. Thorsson et al. 1994). 

For the state this is rigorously true only for isotropic system (see Anderson & 

Morel 1961, Muzikar 1980, and Page 1995). 



2. The spectra of seven pulsars, including the Crab pulsar (PSR 0531+21), can only 
be fitted by a power-law-type spectrum, or by a blackbody spectrum with very 
high effective temperature and effective areas much smaller than the neutron 
star surface. Their X-ray emission is predominated by magnetospheric emission. 
Therefore, the temperatures, determined from the spectral fits, are probably too 
high. Pulsars of this type are marked with dots. 

3. Finally, there are four pulsars, i.e., 0833-45 (Vela), 0656-M4, 0630-M8 (Gemin- 
ga), and 1055-52, allowing two-component spectral fits. The softer blackbody 
component is believed to correspond to the actual surface emission of the neutron 
star, while the harder blackbody (or power law) component may be due to 
magnetospheric emission. These pulsars are marked with squares. 

An overview of the observed luminosities and pulsar ages is contained in Table 5. 



Table 5. Luminosities, L, and spin-down ages, r, of pulsars 



Pulsar 


Name logr [yr] 


logL [erg/s] 


Reference 




not enough data available for spectral analysis 


1706-44 


4.25 


32.8 + 0.7 


Becker et al. (1992) 


1823-13 


4.50 


33.2 + 0.6 


Finley & Ogelman (1993) 


2334-1-61 


4.61 


33.1 + 0.4 


Becker (1993) 




powor-law-type spectra or spectra with only a high temperature component 


0531+21 


Crab 3.09 


33.9 ± 0.2 


Becker & Aschenbach (1995) 


1509-58 


SNR MSH 15-52 3.19 


33.6 + 0.4 


Seward et al. (1983), Trussoni et al. (1990) 


0540-69 


3.22 


36.2 + 0.2 


Finley et al. (1993) 


1951+32 


SNR CTB 80 5.02 


33.8 + 0.5 


Safi-Harb & Ogelman (1995) 


1929+10 


6.49 


28.9 + 0.5 


Yancopoulos et al. (1993), Ogelman (1995) 


0950+08 


7.24 


29.6 + 1.0 


Seward & Wang (1988) 


J0437-47 


8.88 


30.6 + 0.4 


Becker & Triimper (1993) 




spectrum dominated by a 


soft component 


0833-45 


Vela 4.05 


32.9 + 0.2 


Ogelman et al. (1993) 


0656+14 


5.04 


32.6 + 0.3 


Finley et al. (1992) 


0630+18 


Geminga 5.51 


31.8 + 0.4 


Halpern & Ruderman (1993) 


1055-52 


5.73 


33.0 + 0.6 


Ogelman & Finley (1993) 



4-2. Impact of Softening of Pion Exchange Modes and Intermediate Pion Decay on 
Cooling 

In Figs. 7-9 we show the impact of the medium effects on the cooUng of neutron stars 
of different masses. Fig. 7 compares the mass dependence of the neutrino cooling 
rates Lj^/Cy associated with MU-FM79 and MU-VS86 for non-superfluid matter. 
For the soHd curves, the neutrino emissivity in pion-condensed matter (threshold 
density = 3no) is taken into account according to eq. (9) and (1) with the param- 
eter a = \/2. The dashed curves correspond to the case where no pion condensation 
is allowed. As one sees the medium polarization effects included in MU-VS86 may 
result in three order of magnitude increase of the cooling rate for the most mas- 
sive stars. Even for stars of a low mass the cooling rate of MU-VS86 is still few 
times larger than for MU-FM79 because even in this case the more efficient rate is 
given by the reactions shown by the right diagram in Fig.l. The cooling rates of the 
I.SMq mass models with and without pion condensate differ only by a factor of 5. 




0.5 1.0 1.5 2.0 



Fig. 7. Cooling rate due to neutrino emission as a function of star mass for a representative 
temperature of T = 3 x 10* K. The solid curves refer to cooling via MU-FM78+PU and 
MU-VS86+PU in pion-condensed (threshold density ric = 3no) matter. The dashed curves 
refer to cooling without pion condensation. Superfluidity is neglected here. 

One sees that the cooling rates for the FM79 result vary significantly (by about 
three orders of magnitude!) only over a rather small mass range from the models 
without pion condensation to the ones with pion condensation. The cooling rates are 
rather independent of star mass below the critical mass value (i.e., M = 1.63 M©) 
at which the transition into the pion-condensed phase occurs, and show a flattening 



behavior above. The resulting cooHng behavior of neutron stars of several selected 
masses is shown in the left graph of Fig. 8. These curves are quite similar for 
models without pion condensation. Stars which are sufficiently heavy such that a 
pion condensate can develop in their cores evolve along cooling tracks that are 
rather different from the former, provided the stars are older than about 10 years. 
We illustrate this for two representative star masses, M = 1.7 and M = 1.9 M©. 
Depending on star mass, the resulting photon luminosities are basically either too 
high or too low to account for the bulk of observed pulsar luminosities, which tends 
to be a general feature of cooling calculations no matter what kind of enhanced 
cooling mechanism are being studied (cf. Schaab et al. 1996a). 
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Fig. 8. Cooling of non-sup erfluid neutron star models of different masses constructed for 
the HV EOS. The two graphs refer to cooling via MU-FM79+PU (left) and MU-VS86+PU 
(right). In both cases, pion condensation is taken into account for the solid curves where 
n > = 3no. The dashed curves in the right graph refer to the cj^ parametrization 
without pion condensation. The observed luminosities are labeled in Fig. 14. 

The situation changes if the modified Urea process with the medium modifi- 
cations of the pion-exchange interaction, appropriate vertex corrections, and the 



radiation from the intermediate states (MU-VS86) is included. Now the cooHng 
rates vary smoothly with star mass (see Fig. 7) such that the gap between standard 
and enhanced cooling is washed out. More quantitatively, by means of varying the 
star mass between 1.0 and 1.6 Mq, one achieves agreement with a large number of 
observed data points, see right graph of Fig. 8. This is true for both choices of the 
parametrization, independently whether a pion condensation can occur or not. 
The two parametrizations differ only in the range which is covered by the cooling 
curves. The only pulsars which do not agree with the cooling curves are the hot 
PSR's 1055-52 and 1951+32 and the rather cold object PSR 1929+10. The high 
luminosities of the former two may be due to internal heating processes, which leads 
to delayed cooling for star ages t > 10^ yr (see for example Shibazaki & Lamb 1989, 
Sedrakian & Sedrakian 1993, Reisenegger 1995, Van Riper ct al. 1995, Schaab ct al. 
). To achieve agreement with the extremely low luminosity object PSR 1929+10, 
the inclusion of strong magnetic fields in the atmosphere (see Van Riper 1991) or 
other fast cooling processes, like direct Urea etc. seems necessary (cf. Fig. 11 and 
Schaab et al. 1996a). 

Next we shall study the modified Urea process in superfluid matter. Results are 
displayed in Fig. 9. The processes NPBF and PPBF, which will be discussed later, 
have been artificially forbidden. The value of the neutron ^P2 gap from Amundsen & 
0stgaard (1985) is rather large such that the difference between slow (MU-FM79, 
left graph) and intermediate (MU-VS86, right graph) cooling is almost indistin- 
guishable. Superfluidity strongly suppresses the neutrino emission rates from stars 
containing a superfluid core, which thus delays cooling. As it was the case for stan- 
dard cooling, the thermal evolution of superfluid stars is rather insensitive against 
variations of the star's mass. However, as soon as charged pions condense in the core 
of a neutron star, superfluidity disappears and the much higher neutrino emission 
rates cause the cooling curves to drop down by more than two orders of magnitude. 
For the model which allows for pion condensation, this occurs when the star's mass 
is varied from 1.6 to 1.7 Mq. In principle, it appears possible to achieve agreement 
between the body of observed cooling data and the cooling simulations by properly 
varying the star's mass in an extremely narrow range. Surely, such a procedure 
is not very satisfying because it requires an extreme fine-tuning of the mass of a 
star. As we shall see below, the superfluid state allows for other neutrino-emitting 
processes which decisively influence the cooling of stars and thus provide other, less 
stringent conditions by means of which agreement with the observed data can be 
achieved. 




4-3. Superfluid Pair Breaking and Formation 

To demonstrate the relevance of the superfluid pair breaking and formation pro- 
cesses for the cooling of neutron stars, in the first step, we implement these two 
processes into standard cooling calculations. In these calculations the FM79 result 
for the modified Urea process is used. No further processes, like DU or PU is taken 
into account here. The results are shown in Fig. 10. One sees that NPBF and PPBF 
reduce the photon luminosity of neutron stars significantly for star ages between 
about 1 < t < 10^ yr. The acceleration of the cooling is mostly due to the super- 
fluid transition of the neutrons in the crust since the gap energy for this transition 
is 1.13 MeV compared to 0.25 MeV for the protons in the core (see Tab. 4) and 
the value of A/T at T = Tc is 0.57 compared to 0.12 for the neutrons in the core 
(see Sec. 2.2). One can also see that the importance of the NPBF process increases 
with decreasing gap energy of the ^ So-neutron-pairing. For t > 10^ yr the influence 
of NPBF and PPBF becomes rather weak and cooling via photon emission from 
the star's surface becomes the dominant process. As a result, the solid and dashed 
curves in Fig. 10 are almost identical for t > 10^ yr. 




Fig. 10. Cooling of neutron stars of mass M = 1.4Mq constructed for the HV EOS. 
The curves correspond to standard cooling (MU-VS86 and PU ignored) with (solid curve) 
and without (dashed curve) NPBF+PPBF. The gap energy of the ^ So-neutron-pairing is 
arbitrarily enhanced (dotted curve) or reduced (dashed-dotted curve), respectively by a 
factor of 2 in order to demonstrate the dependency of the rates on the gap energy. The 
observed data are labeled in Fig. 14. 



4-4- Final Results and Influence of Equation of State 

We turn now to cooling simulations where the MU, NPBF, PPBF, DU and PU take 
place simultaneously. Fig. 11 shows the cooling tracks of stars of different masses, 
computed for the HV equation of state. The DU process is taken into account in the 
right graph, whereas it is neglected in the left graph. The solid curves refer again 
to the cD^-parameterization with phase transition to a pion condensate, the dashed 
curves to the one without phase transition. For masses in the range between 1.0 
and 1.6 Mq, the cooling curves pass through most of the data points. We again 
recognize a photon luminosity drop by more than two orders of magnitude of the 
1.7 Mq mass star with pion condensate, due to vanishing of the superfluidity as 



a consequence of pion condensation. This drop is even larger if the DU is taken 
into account (right graph, this allows to account for the photon luminosity of PSR 
1929+10). 

A comparison with the observed luminosities shows that one gets quite good 
agreement between theory and observation if one assumes that the masses of some 
of the miderlying pulsars are different from the canonical value, M = 1.4 Af©. 
Since the masses of these pulsars are not known, no further conclusions about the 
actually operating cooling mechanism of these pulsars can be drawn yet. Future 
mass determinations of these objects will change the situation. The only two stars 
whose photon huninosities cannot be accounted for by the new processes studied 
here are the rather hot pulsars PSR 1055-52 and PSR 1951+32. Again, as mentioned 
above, their temperatures may be explained by internal heating. 
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Fig. 11. Cooling of neutron stars with different masses constructed for the HV EOS. The 
cooling processes are MU-VS86, PU (only solid curves), NPBF, PPBF, and DU (only in 
the right graph). The dashed curves refer to the M = 1.7 and 1.9Mq models without pion 
condensate. The observed data are labeled in Fig. 14. 



To explore the impact of somewhat smaller superfluidity gaps on cooling, we 
arbitrarily reduced the gap by a factor of 6. Its major effect is a drop in the cooling 
curves of the 1.0 to 1.6 Mq mass stars and of the 1.7 and 1.9 Mq mass stars without 
pion condensate (dashed curves). This would lead to better agreement with the low- 
temperature pulsar 0630+18 (Geminga) (see Fig. 12). 




Fig. 12. Same as left graph of Fig. 11 but with the original '^P2 gap by Amundsen & 
0stgaard (1985) reduced by a factor of 6. The observed data arc labeled in Fig. 14. 



So far we have selected a particular model for the equation of state of neutron 
star matter, i.e. HV, and studied the dependence of the cooling curves on the star's 
mass. The influence of different models for the equation of state on cooling is shown 
in Figs. 13-15. Fig. 13 displays the cooling behavior computed for a representative 
non-relativistic model for the equation of state, UV14+UVII. The general qualita- 
tive features obtained for HV carry over to this case, except that the temperature 
drop occurs now at a lower mass threshold. For that reason the 1.1 Mq mass star 
is already anomalously cold. Because of the relatively high neutron fraction of the 



UVi4+UVII compared to relativistic EOS's the superfluid phase does not reach to 
the center of the 2.\Mq mass star in this particular model of Amundsen & 0stgaard 
(1985). This causes the luminosity drop to occur also for the model without pion 
condensate (dashed curve) . The high neutron fraction, or identically the low proton 
fraction, has also the effect that the DU process cannot occur. 
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Fig. 13. Cooling of neutron stars of different masses constructed for the UV14-I-UVII 
EOS. The cooling processes are MU-VS86, PU (only solid curves), NPBF, and PPBF. 
The observed data are labeled in Fig. 14. 

The cooling curves constructed for RBHF, shown in Fig. 14, are rather similar to 
those of UV14+UVII though the underlying microscopic theories employed to deter- 
mine these two models for the equation of state are completely different. We recall 
that the former, RBHF, is based on a parameter-free relativistic Brueckner-Hartree- 
Fock calculation while the latter, UV14+UVII, was derived in the framework of a 
non-relativistic many-body variational calculation. In both cases, however, dynami- 
cal nucleon-nucleon correlations are taken into account, which - as a general feature 
- tend to soften the equation of state at intermediate nuclear densities, both in the 



relativistic as well as in the non-relativistic case (Fig. 4). For that reason both of 
these models for the equation of state lead to similar slopes of star mass as a func- 
tion of central star density (see Fig. 5) and similar cooling behaviors. With respect 
to the differences between RBHF and UV14-I-UVII, the former is harder than the 
latter. Therefore the magnitude of the temperature drop for increasing star masses 
is smaller. For both of these models, pion condensation set in at rather low neu- 
tron star masses, causing already the rather light neutron stars (M ~ 1.1 Mq) to 
cool very rapidly. In contrast to the UV14-I-UVII EOS, the RBHF has large proton 
fractions at high densities causing the superfluid phase to reach to the center of the 
star for all star masses and allowing for the DU process. The cooling tracks for the 
(1)^ -parameterization without pion condensation and the ones with the DU process 
arc similar to the cooling tracks for the HV EOS. Thus, we can restrict ourselves 
to the pion condensation case without DU process. 
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Fig. 14. Cooling of neutron stars of different masses constructed for the RBHF EOS. The 
cooling processes are MU-VS86, PU, NPBF, and PPBF. The observed data are labeled in 
Fig. 14. 



Besides HV, we employ a second model, labeled G300, for the equatfon of state 
derived in the framework of relativistic mean field theory. Both equations of state 
are rather similar, as can be seen in Figs. 4 and 5, and so are the corresponding 
cooling curves. Fig. 15. Again only the pion condensation case without DU process 
is studied. 
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Fig. 15. Same as Fig. 14 but for the G^°° EOS. The observed data are labeled in Fig. 14. 



5. Summary 

Two recalculated neutrino-emitting processes in dense neutron star matter, i.e., the 
modified Urea process and the superfiuid pair breaking and formation processes, 
have been implemented in numerical cooling simulations of neutron stars. Our major 
finding is that standard cooling processes supplemented with the medium effects 
enables one to achieve theoretical agreement with a large fraction of the observed 
pulsar luminosities, which — subject to uncertainties in the equation of state and 



transport properties of superdense matter — seems to be a problem for the standard 
cooling scenario alone (of. Schaab ct al. 1996a). In particular it is possible to account 
for the low temperatures of PSRs 0833-45 (Vela), 0656+14, and 0630+18 (Geminga). 
Secondly, the new processes studied here provide a crossover between standard 
cooling and the so-called fast cooling scenarios of neutron stars and thus can be 
viewed as "intermediate" cooling scenarios. 
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